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Assuming the Majorana nature of massive neutrinos, we generalize the Friedberg-Lee 
neutrino mass model to include CP violation in the neutrino mass matrix M v . The most 
general case with all the free parameters of M v being complex is discussed. We show 
that a favorable neutrino mixing pattern (with 9 12 ~ 35.3°, 9 23 = 45°, 9 13 ^ 0° and 
S = 90°) can naturally be derived from M u , if it has an approximate or softly- broken fi-r 
symmetry. We also point out a different way to obtain the nearly tri-bimaximal neutrino 
mixing pattern with <5 = 0° and non-vanishing Majorana phases. 
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Recently, a novel neutrino mass model has been proposed by Friedberg and Lee 
(FL) P The neutrino mass operator in the FL model is simply given by 

^-mass =a{pT~ V a) { V r ~ V y) + b ~ "e) ( V a ~ V e) + C { V e ~ V r) Oe ~ V r) 

+ m (Ve v e + + V r V T ) , (1) 

where the parameters a, b, c and m are all assumed to be real, and the charged- 
lcpton mass matrix is taken to be diagonal. A salient feature of £„_ mass is its partial 
gauge-like symmetry; i.e., its a, b and c terms are invariant under the transformation 
v a — » v a +z (for a = e, /i, t) with z being a space-time independent constant element 
of the Grassmann algebra.^ From Eq. (1), one can directly write down the neutrino 
mass matrix: 



M v = m 1 + \ -b a + b -a . (2) 





Two interesting features can be inferred from the diagonalization of M v . First, the 
neutrino mass matrix takes a magic form,2 in which the sums of rows and columns 
are all equal to m . The unitary matrix used to diagonalize M v must have one 
eigenvector with three equal components \j\f2> . Second, when b = c holds, it is 
very easy to check that the neutrino mass operator £j,_ mass has the exact /z-r 
symmetry (i.e., £„_ mass is invariant under the exchange of /i and r indices).^ In 
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addition, one may consider to remove one degree of freedom from £„_ mass or M v 
(for instance, by setting c = 0)P 

To include CP or T violation into the FL model, one may insert the phase 
factors e ±4 '' into Eq. (1) by replacing the term c (77 e — V T ) [y e — v T ) with the term 
c (e~ l71 v e — v T ) (e + '"V e — ^ T )- The resultant neutrino mass matrix is no longer 
symmetric, hence it describes Dirac neutrinos instead of Majorana neutrinos. How- 
ever, in most of the realistic models, the Majorana nature is preferable to the Dirac 
nature of neutrinos. Hence, in this work, we aim to generalize the FL model to 
include CP and T violation for massive Majorana neutrinos. 

Let us start from the generic analysis with all the parameters of M v in Eq. (2) 
being complex. For Majorana neutrinos, M v is symmetric and can be diagonalized 
by the transformation V' M V V* — Diag{rn l5 m 2 , m 3 }, in which m i (for i — 1,2,3) 
stand for the neutrino masses. After a straightforward calculation, the neutrino 
mixing matrix V turns out to be 



V 



— = cos - 

V6 2 



— = cos = sm — e 

V6 2 y/2 2 



;= cos - 

V V6 2 



—= sm — e 

y/2 2 
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-i6 



— = sm — e 

y/6 2 



— == cos = sm — e 
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-iS 



= cos = sm — e 

y/2 2 V6 2 



-id 



(e tp 0\ 

e la 
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where the explicit expressions of 9 and S are 

s/X^ 



tan# = 



F 2 



r X 

tano = — 



(3) 



(4) 



and the definitions of X, Y and Z can be found in Ref. HI Furthermore, three mass 
eigenvalues of M v and two Majorana phases of V are found to be 



r l cos 2 --T 2 sin6'e 



-id 
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and 
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Again, the expressions of have also been listed in Ref. |H 
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We proceed to consider two special but interesting scenarios of the generalized 
FL model and explore their respective consequences on three neutrino mixing angles 
and three CP-violating phases. 

Scenario (A): a and m are real, and b — c* are complex. Note that the 
fi-r symmetry of M v is softly broken in this case, because |6| = |c| holds. By 
using the generic results given in Eqs. (3)- (6), one can easily arrive at tan# = 
v / 3Im(6)/[™ +a + 2Rc(6)], S = 90°, 



m 1 = y [m + a + 2Re (6)] 2 + 3 [Im (6)] 2 — a + Re (6) , m 2 = m , 

m 3 = \J [rn + a + 2Re {b)f + 3 [Im (b)f + a - Re (b) , (7) 

together with p = a = 0. Comparing our results with the well-known standard 
parametrizationJSlwe immediately obtain sin0 12 = 1/ (V2 + cos6»), sin(9 23 = l/y/2 , 
and sin6* 13 = 2/ k/6sin(0/2)l . The leptonic Jarlskog parameter J, which is 
a rephasing-invariant measure of CP violation in neutrino oscillations,® reads 
J = sin6'/(6v / 3)- If = 0° holds, the tri-bimaximal neutrino mixing pattern (with 
tan(9 12 = l/y/2 or 12 » 35.3°, 6 23 = 45° and 6 13 = 0°) 7 will be reproduced. One 
can see that the soft breaking of \i-t symmetry leads to both 6 13 ^ 0° and J ^ 0, 
but it does not affect the favorable result 9 23 — 45° given by the tri-bimaximal 
mixing pattern. On the other hand, sin# 12 w l/v3 is an excellent approximation, 
since 8 must be small to maintain the smallness of 9 13 . In view of 9 13 < 10°f^we 
obtain 9 < 24.6° and J < 0.04. It is possible to measure J - O(10" 2 ) in the 
future long-baseline neutrino oscillation experiments. The neutrino masses in sce- 
nario (A) rely on four real model parameters m , a, Re(6) and lm(b). Thus it is easy 
to fit the neutrino mass-squared differences Am 21 = (7.2 . . . 8.9) x 10~ 5 eV 2 and 
Am§ 2 = ±(2.1 . . .3.1) x 10 -3 eV 2 P Such a fit should not involve any fine-tuning, 
because (a) the number of free parameters is larger than the number of constraint 
conditions and (b) three neutrino masses have very weak correlation with three mix- 
ing angles. A detailed numerical analysis can be found in Ref. [H and a remarkable 
feature is that only the normal neutrino mass hierarchy (m 1 < m 2 < m 3 ) is allowed 
in this scenario. 

Scenario (B): a, b and c are all real, but m is complex. By using Eqs. (3)-(6), 
we obtain tan 8 — y/3 (b — c) / (b + c — 2a) , and 



m 1 = y [m_ + Re (m )] + [Im (m )] 2 , m 2 = \m \ , 

m 3 = J [m + + Re (m )] 2 + [Im (m )] 2 , (8) 



where m_(_ = (a + b + c) ± V a 2 + b 2 + c 2 — ab — ac — be . Two Majorana phases p 
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and a are given by 



tan 2p 



(m + — m_) Im (m ) 



tan 2(7 



nip + m + m_ + (m + + m_) Re (m ) ' 
m + Im (m ) 



(9) 



TOq + ?7i + Re (m ) 



Although p and er have nothing to do with the behaviors of neutrino oscillations, 
they may significantly affect the neutrinoless double-beta decayP^ Comparing our 
formulae with the standard parametrization, we arrive at sin# 12 = l/v2 + cos 9, 
sin6> 23 = \Jl + cos 9 - V3 sin 9/^/2 (2 + cos6>), sin6» 13 = 2/V6 sin(0/2) together 
with 6 — 0° for the Dirac phase of CP violation. The results for 9 12 and 6* 13 in this 
scenario are the same as those obtained in scenario (A) , but the Jarlskog parameter 
J is now vanishing. Because of the p-r symmetry breaking, 9 23 may somehow devi- 
ate from the favorable value # 23 — 45°. Given 9 < 24.6° corresponding to 9 13 < 10°, 
# 23 is allowed to vary in the range 37.8° < 9 23 < 45°. The neutrino masses depend 
on five real model parameters a, b, c, Re(m ) and Im(m ). Hence there is sufficient 
freedom to fit two observed neutrino mass-squared differences Am^ and Am^. 
Our careful numerical analysis, which has been done in Ref. [H shows that both 
normal (m 1 < m 2 < m 3 ) and inverted (m 3 < m 1 < m 2 ) neutrino mass hierarchies 
are allowed, and the Majorana phases (p, a) are less restricted in scenario (B). 

Although our discussions about the generalized FL model are restricted to low- 
energy scales, it can certainly be extended to a superhigh-energy scale (e.g., the 
GUT scale or the seesaw scale). In this case, one should take into account the 
radiative corrections to both neutrino masses and fl avor mixing parameters when 
they run from the high scale to the electroweak scaleff^ 

We conclude that the p-r symmetry and its slight breaking are useful and sug- 
gestive for model building. We expect that a stringent test of the generalized FL 
model, in particular its two simple and instructive scenarios, can be achieved in 
the near future from the neutrino oscillation and neutrinoless double-beta decay 
experiments. 
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